arXiv: 1507.04739v2 [math.CO] 5 Nov 2015 


Counting matchings in irregular bipartite graphs and random lifts 


M. Lelarge* 


Abstract 

We give a sharp lower bound on the number of matchings of a given size in a bipartite 
graph. When specialized to regular bipartite graphs, our results imply Friedland’s Lower 
Matching Conjecture and Schrijver’s theorem proven by Gurvits and Csikvari. Indeed, our 
work extends the recent work of Csikvari done for regular and bi-regular bipartite graphs. 
Moreover, our lower bounds are order optimal as they are attained for a sequence of 2-lifts of 
the original graph as well as for random n-lifts of the original graph when n tends to infinity. 

We then extend our results to permanents and subpermanents sums. For permanents, 
we are able to recover the lower bound of Schrijver recently proved by Gurvits using stable 
polynomials. Our proof is algorithmic and borrows ideas from the theory of local weak 
convergence of graphs, statistical physics and covers of graphs. We provide new lower bounds 
for subpermanents sums and obtain new results on the number of matching in random ?r-lifts 
with some implications for the matching measure and the spectral measure of random n-lifts 
as well as for the spectral measure of infinite trees. 


1 Introduction 


Recall that a n x n matrix A is called doubly stochastic if it is nonnegative entrywise and each 
of its columns and rows sums to one. Also the permanent of a n x n matrix A is defined as 

n 

per(A) = X 

ffSiSn *=1 

where the summation extends over all permutation a of {1,... ,n}. The main result proved in 
m is the following theorem: 

Theorem 1. (Schrijver For any doubly stochastic n x n matrix A = we define 

A = (fiij = aij( 1 — ajj)) and we have 

per(A) > JJ(1 - a itj ). (1) 

It is proved in ED 122] that this theorem implies: 

Theorem 2. Let A be a non-negative n x n matrix. Then, we have 

In peg A) > max (1 — x % j) ln(l — x t ,) + X{, In ( -LL ] , (2) 

n . . \ Xi j ] 

l,J x ' 
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with the convention In ^ = 1 and where M„ iti is the set of n x n doubly stochastic matrices. 


Clearly applying Theorem [2] to A with = aij, we get Theorem [T| back, so that both 
theorems are equivalent. In [T9 , 20] , L. Gurvits provided a new proof of these theorems using 
stable polynomials, see also [25]. Our main new result is a generalization of Theorem [2] to 
subpermanent sums, see Theorem [4] below. Our proof is very different from those derived in 
and borrows ideas from the recent work of Csikvari [14] for matchings in regular 
bipartite graphs. Permanents and subpermanent sums can be interpreted as weighted sums 
of matchings in complete bipartite graphs and the extension of the approach in [14] to this 
general framework is the main technical contribution of our work. Interestingly, the obstacles 
to overcome are computational in nature and we present a very algorithmic solution. As a 
byproduct, we also obtain new results on the number of matchings in random lifts, see Theorem 
[5]below. We compute the limit of the matching generating function (i.e. the partition function of 
the monomer-dimer model) for a sequence of random n-lifts of a graph G as an explicit function 
of the original graph G. This result has also some implications for the spectral measure and the 
matching measure of random lifts, see Theorems [7] and [8] 

As a consequence of Theorem [TJ Schrijver shows in [37J that any d-regular bipartite graph 
with 2n vertices has at least 

{d- l)^ 1 ! 
d d ~ 2 ) 


perfect matchings (a perfect matching is a set of disjoint edges covering all vertices). For each 
d, the base (d— 1 ) d ~ l /d d ~ 2 in ([3]) is best possible [JTj. Similarly, our Theorem [4] shows that the 
right-hand term in ([2]) is best possible in the following sense: for any n x n matrix A, we show 
that there exists a sequence of growing matrices Ai obtained by taking successive 2-lifts of A 
(see Definition © such that its permanent grows exponentially with its size at a rate given by 
the right-hand term in ([2]) . We refer to Theorem [7] for a precise statement and similar results 
for subpermanent sums. 

In [16] . Friedland, Krop and Markstrom conjectured a possible generalization of ([3]) which 
is known as Friedland’s lower matching conjecture: for G a d-regular bipartite graph with 2 n 
vertices, let mi~{G) denote the number of matchings of size k (see Section 12.11 for a precise 
definition), then 


where p = An asymptotic version of this conjecture was proved using Theorem Q] in [21,12 2], 
A slightly stronger statement of the conjecture was proved by Csikvari in [I7j and we extend it 
to cover irregular bipartite graphs, see Theorem [3] 

We state our main results in the next section. In Section [3] we summarize the main ideas of 
the proof and describe related works. We also give some implications of our work for extremal 
graph theory and for the spectral measure and matching measure of random lifts. Section [4] 
contains the technical proof. We first summarize the statistical physics results for the monomer 
dimer model in Section 1441 Then, we study local recursions associated to this model in Section 
14.21 The results in this section build mainly on previous work of the author [28]. In Section 14.31 
we show how an idea of Csikvari m using 2-lift extends to our framework and connect it to 
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the framework of local weak convergence in Section 14.41 We use probabilistic bounds on the 
coefficients of polynomials with only real zeros to finish the proof in Section 14.51 In Section 14.61 
we provide the details needed for random lifts. 


2 Main results 

We present our main results in this section. The results concerning lower bounds for the number 
of matchings given in Section 12.11 are implied by those in Section 12.21 concerning lower bounds 
for permanents. The reader only interested in the most general result concerning lower bounds 
might jump directly to Section T2.2I where Theorem [4] is stated in a self-contained manner. Note 
that the proof will use notations introduced in Section 12.11 Section 12.31 contains our results 
on random lifts. They are independent from the other results but require some notations from 
Section 12.11 


2.1 Lower bounds for number of matchings of a given size 

We consider a connected multigraph G = (V,E). We denote by v(G) the cardinality of V: 
v(G) = \V\. We denote by the same symbol dv the set of neighbors of node v € V and the set 
of edges incident to v. A matching is encoded by a binary vector, called its incidence vector, 
B = (B e , e £ E) € {0, 1} E defined by B e = 1 if and only if the edge e belongs to the matching. 
We have for all v € V, Yle&dv P e < 1. The size of the matching is given by Yh e p e- We will also 
use the following notation e € B to mean that B e = 1, i.e. that the edge e is in the matching. 
For a finite graph G, we define the matching number of G as v(G) = max{^ e B e } where the 
maximum is taken over matchings of G. 

The matching polytope M(G) of a graph G is defined as the convex hull of incidence vectors 
of matchings in G. We define the fractional matching polytope as 


FM(G) 



(5) 


We also define the fractional matching number v*{G) = max xgFM ( G ) Yh e x e — It is 

well-known that: M(G ) = FM(G ) if and only if G is bipartite and in this case, we have 
i/(G) = i/*(G). 

For a given graph G, we denote by rrik(G ) the number of matchings of size k in G (mo(G) = 
1). For a parameter z > 0, we define the matching generating function: 


u{G) 

p g(z) = m k(G)z k . 

k =0 

In statistical physics, the function lnPcOz) is called the partition function. We define by Mfc(G) 
the convex hull of incidence vectors of matchings in G of size k and similarly for 0 < t < v*(G): 


FM t (G) 


X £ 


x e > 0, 


E 

e£dv 


X e < 



( 6 ) 
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If G is bipartite, we have M k (G) = FM k [G). Theorem [3] below deals with bipartite graphs but 
Theorem 0 deals with general graphs and fractional matching polytopes will be needed there. 

For any finite graph G , we define the function Sq : FM(G) —» R by: 

Sc( x ) = Y -x e In x e + (l- x e ) ln(l - x e ) - Y ( 1 - Y x e j In f 1 - Y x e j . (7) 

e€E v€V \ e£dv / \ e£dv / 

The function Sq is concave on FM(G) by Proposition 1131 

Definition 1. Let G be a graph with no loops. Then H is a 2-lift of G ifV(H) = V(G) x {0,1} 
and for every (u,v) G E(G), exactly one of the following two pairs are edges of Ft : ((it, 0), (v,0)) 
and (( u , 1), (v, 1)) G E(H ) or ((it, 0), (v, 1)) and ((it, 1), ( v , 0)) G E(H). If (it, v) E(G), then 
none of ((it, 0), (v, 0)),((it, 1), (u, 1)), ((it, 0), (u, 1)) and ((it, 1), (u, 0)) are edges in H. 


Note that a bipartite graph G has no loops. 

Theorem 3. For any finite bipartite graph G, we have for z > 0, 


In Pa(z) > max 
xeM(G) 


lnz + 5|(x)| 


( 8 ) 


We /lane for all k < v(G), 


m k (G) > b v i G) k {k/v{G)) exp ( max S§{x) ) , 

Vxe M k {G) ) 

where b n k (p) is the probability for a binomial random variable Bin(n,p) to take the value k, 
i.e. b ntk (p) = (k)p fc ( 1 — p) n ~ k . Moreover, there exists a sequence of bipartite graphs {G n = 
(V n , E n )} n£ ?$ such that Go = G, G n is a 2-lift of G n -1 for n > 1 and for all z > 0, 


lim } In P Gn (z) = —max 
n-> °o v(G n ) v{G ) xSM(G) 


Y X ^j lnz + S G ( x ) | 


Consider the particular case where G is a d-regular bipartite graph on 2n vertices. In this 
case, we have v(G) = n and we can take x* = — for all e G E so that x* G M k (G) and we have 


S8(*-*)=n 


(p In (p) + ( d ~ ln ( X ~~ d) ~~ 2< ' 1 ~~ ^ ln ^ “ ’ 


with p = y . We see that we recover the first statement in Theorem 1.5 of [33]. In particular, for 
k = n, i.e. p = 1, we recover Q and for k < n, as explained in [13], we slightly improve upon 
(HD- Note that in this particular case, we have m n {G ) < ( d\) n / d by a result of Bregman [12] (see 
also m for upper bounds for m k (G) with k <n). 

Taking z = 1 in (|8l) , we obtain the following bound on the total number of matchings: 
Proposition 1. For any bipartite graph G, we have: 

PG) 

m k (G) > exp ( max S8(x) ) . 

t'o P Vx£M(G) V 
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2.2 Lower bounds for permanents 

In this section, we extend previous results to weighted graphs. We state our results in term of 
permanents. Let A be a non-negative n X n matrix. We denote by M n the set of such matrices. 
Recall that the permanent of A € M n is defined by 

n 

per(d) = y n^M- 
( 7GS n 2=1 

We define by M n n the set of n x n doubly stochastic matrices: 


M n n = l A, 0 < a, 


i,3) 


E- 


E- 


= 1 > C M n 


For 1 < k < n, let per fc (d) be the sum of permanents of all k x k minors in A and per 0 (d) = 1. 
per fc (d) is called the k -th subpermanent sum of A. We define M n ^ the set of n x n non-negative 
sub-stochastic matrices with entrywise Li-norm k: 


M U) k = 


a, o ^ *y ^ * 


H,j < 


E' 

j 


H,j < 


i.E' 


J i,3 


= k > C M n 


h3 


We also define the set of substochastic matrices: 


Mn,< — \ A, 0 < O’iji 'y ^ Q'ij — 1) 'y ^ Qjjj ^ 1 ^ C M n . 

{ i 3 

We define the function S B : M n x M n .< ->1U {—oo} by 
S b (A,x) = y Xjj In — + (1 - Xjj) ln(l - Xjj) 


(9) 


h3 


^ ~ J2 Xi ’^j ln ^ ^ - 5Z Xi ^ ln ^ 


with the convention ln jj = 1. First note that if A is the incidence matrix of a bipartite graph 
G , and x is such that there exists Ojj = 0 and Xij > 0, then S B (A,x ) = —oo. Moreover 
if x has only non-negative components corresponding to edges of the graph G, then we have 
S B (A,x) = S B (x) as defined in fl7|) with a slight abuse of notation: the zero components (on 
no-edges of G) of x as argument of S B (A. x) are removed in the argument of S B (x). Note that 
x i —> S b (A,x.) is concave on M n ^< (since S B is concave on FM(G) by Proposition fl3l) . 

Definition 2. Let A be a non-negative n x n matrix. Then B is a 2-lift of A if B is a 2n x 2 n 
non-negative matrix such that for all i,j E {1,.. . n}, either bij = bi +n j +n = a h] and bij +n = 
b ,. n.j — 0 or bij+n — b, - ii.j — Oj.y and b^^j — ^-j-nj+n — 0* 
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Theorem 4. Let A be a non-negative n x n matrix. Let v(A) = max{fc, per fc (A) > 0}. For all 
k < v{A), we have 


per fc (A) > b v ( A ) tk (k/v(A))exp 


I max 

\x.£M n ,k 



(ii) 


where b n ^{p) = Q!')p fc (l — p) n k . Moreover, there exists a sequence of matrices {A^ € M 2 e n }e^ 
such that Ao = A, Ai is a 2-lift of Ae_\ for £ > 1 and for all z > 0, 


lim —t In 

e->oo 2 £ 


G(A e ) 

P er k( A l) 

k =0 


Z = 


max 

x6M„,< 



In z + S B (A,x 


If k = n in ED, we recover Theorem [2] which is equivalent to Theorem [T] Note that if 
v(A) < n, then per(A) = 0 and the lower bound in Theorem [2] is equal to —oo. Indeed if 
per(A) = 0, then if x € M n ^ n is a permutation matrix then there exists i,j such that a,;j = 0 

and Xij > 0 so that In ^^4^ = —oo. The claim then follows from the Birkhoff-von Neumann 
Theorem which implies that any doubly stochastic matrix can be written as a convex combination 
of permutation matrices. Also, results presented in Section [2.11 follow by taking for the matrix 
A, the incidence matrix of the graph G. 


2.3 Number of matchings in random lifts 

As in previous section, G = ( V , E ) is a fixed connected multigraph with no loops. A random 
n-lift of G is a random graph on vertex set V\ U V 2 U ■ ■ ■ U K,(g) , where each Vi is a set of n 
vertices and these sets are pairwise disjoint, obtained by placing a uniformly chosen random 
perfect matching between V t and Vj, independently for each edge e = ij of G. We denote the 
resulting graph L n (G). 

Our main result shows that the lower bounds derived in Section 12.11 are indeed attained by 
a sequence of random lifts. More precisely, we have 

Theorem 5. For any finite graph G, we have 

lim -v(L n (G)) 

n—>• 00 Ti 

Vz> 0, lim -\nPj {G) (z) 

where S B (x) is defined in 0. We denote by FPM(G) = FM v * {g) {G) = {x, X eea „ x e = 1} 
the fractional perfect matching polytope of G, then we have 

lim sup -\nm u ( Lji (Q)){L n (G)) < sup S G {x). 
n—>00 n x£FPM(G) 

If, in addition, G is bipartite, then the fractional perfect matching polytope is simply the perfect 
matching polytope PM (G) of G and 

lim -In m v i L , G ML n (G)) = sup S G (x). 
n^-oo n xePM(G) 


= v*{G) a.s. 


max 

xeFM(G) 



In z + S G (x) ^ a.s. 
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In |31| . Linial and Rozenman studied the existence of a perfect matching in L n (G). Note 
that if the number of vertices in G is odd, then in order to have a perfect matching in a n-lifts 
of G, we need to have n even. For n even, they described a large class of graphs G for which 
L n {G) contains a perfect matching asymptotically almost surely. This class contains all regular 
graphs and, in turn, is contained in the class of graphs having a fractional perfect matching, i.e. 
graphs G such that 2v*(G) = v(G). Our result shows that in this case, L n (G ) will contain an 
almost perfect matching (possibly missing o(n) vertices) almost surely. If in addition the graph 
G is bipartite, the number of such matchings is exponential in n. 

In p 8 j, the number of perfect matchings in L n (G), denoted by pm(L n (G)), is studied in the 
limit n —>• oo (along a subsequence ensuring the existence of a perfect mating), where G is a 
graph with a fractional perfect matching. Using the small subgraph conditioning method, an 
asymptotic formula for E[pm(L n (G’))] is computed for any connected regular multigraph G with 
degree at least three. Partial results are also given for E[pm(L n (G)) 2 ] with an explicit formula 
based on a conjecture (proved only for 3-regular graphs). 

Note that we always have v(L n {G)) > nv(G) (by lifting a maximum matching of G). In 
particular, if G has a perfect matching then i/(L n (G)) = nv{G)/ 2 . Hence, Theorem [5] implies 
that for any d-regular bipartite graph G, we have 


lira — lnpm(L n (G)) 
n—> oo Tl 



(d-l) d ~ 1 \ 
d d ~ 2 ) ' 


This result is consistent with [18j . Indeed by Jensen inequality, we always have 


E [In pm(L n (G))] < InE [pm(L n (G))] , 


and our result shows that in the large n limit, the two quantities are asymptotically equal. 
Note that the fact that E [pm(L n (G)) 2 ] ~ E [pm(L n (£?))] 2 to leading exponential order is not 
sufficient to prove this asymptotic equality. 

A similar result to Theorem [5] was shown for permanent in [391140] . It is possible to define 
a random fc-lift for a matrix A (by taking the weighted biadjacency matrix of the random fc-lift 
of the weighted graph associated to the biadjacency matrix A). It is shown in [39, [JOj that 
limsupfc^.^ iE[lnper(L fe (A))] = max ieM „, B “ x i,j) M 1 ~ x i,j) + x i,j ln ( 777 ) • Since [333 

recently showed that for Ak any /c-lift of A, we have per(Afc) < per(A) fc , they obtained a new 
proof of © . 


3 Main ideas of the proof and more related works 

Recall that ln Pg{z) = In Y!k m k{G)z k is called the partition function. A crucial observation 
made by Csikvari [14| is that for any 2 -lift H of a bipartite graph G, we have ln Pg(z) > 

hi Ph(z), see Proposition [ 8 ] in the sequel. Note that a 2-lift of a d-regular graph is still 
d-regular. Starting form any d-regular bipartite graph G, Csikvari builds a sequence of 2-lifts 
with increasing girth. Then Csikvari uses the framework of local weak convergence in order to 
define a limiting graph for the sequence of 2-lifts. In this particular case, the limit is the infinite 
d-regular tree. Using the connection between the local weak convergence and the matching 
measure of a graph developed in QQ|, Csikvari computes a limiting partition function associated 
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to this infinite d-regular tree which is obtained from the Kesten-McKay measure. This limiting 
partition function is then a lower bound for the original partition function In Pg(z) and the 
lower bounds (ED and ED are then easily obtained by properly choosing the parameter z as a 
function of k the size of the matchings that we need to count. 

Our work extends the recent work of Csikvari done for regular and bi-regular bipartite graphs 
to irregular bipartite graphs. We are still building a sequence of 2-lifts with increasing girth and 
still using the framework of local weak convergence. The limiting object is now the universal 
cover of the initial graph G, i.e. the tree of non-backtracking walks also called the computation 
tree (see Section roi for a precise definition). A direct computation of the matching measure 
of this (possibly infinite) tree seems tedious. Here, we depart significantly from the analysis of 
Csikvari. Our approach for the computation of the limiting partition function is based on an 
alternative (more algorithmic) characterization first developed in [1 lj based on local recursions 
on the universal cover of G. In order to express the computations done on the universal cover as 
a function of the original graph G, we rely on results proved by the author in [28] where the local 
recursions are studied on any finite graph. The solution of these local recursions on the universal 
cover (hence a possibly infinite tree) is in correspondence with the solution of the local recursions 
on the initial graph. Since this solution is given by the maximum of a certain “entropy-like” 
concave function defined by 0 on the fractional matching polytope of the original graph, we 
obtain an explicit formula for the limiting partition function and the lower bound in Theorem 
[3] follows. Our approach is then generalized to weighted bipartite graphs in order to get our 
results for (sub-)permanents (Theorem 0 and to random lifts in order to get Theorem [5] 

As explained in the sequel of this section, all the basic ideas used in our proofs were present 
in a form or another in the literature: using lifts for extremal graph theory was one the main 
motivation for their introduction in a series of papers by Amit, Linial, Matousek, Rozenman 
and Bilu m eh Eii; the matching measure and the local recursions already appeared in the 
seminal work of Heilnrann and Lieb [23]; the function Sq defined in 0 is known in statistical 
physics as the Bethe entropy [53]. The main contribution of this paper is a conceptual message 
showing how known techniques from interdisciplinary areas can lead to new applications in 
theoretical computer science. In the next subsections, we will try to relate our results to the 
existing literature and give credit to the many authors who inspired our work. 

3.1 Covers, extremal graph theory and message passing algorithms 

The idea to use lifts to build graphs with extremal properties is not new. In [8], Bilu and Linial 
study 2-lift of d-regular graphs in order to construct infinite families of expanders. They showed 
that the eigenvalues of a 2-lift are the union of the eigenvalues of the original graph and those of 
the signing associated to the 2-lift. They conjectured that every d-regular graph has a signing 
with spectral radius at most 2 \Jd — 1. This conjecture was proved by Marcus, Spielman and 
Srivastava in [32j where they construct bipartite Ramnujan graphs of all degree. 

We can informally state our Theorem [3] as an extremal graph theoretic result: among all 
bipartite graphs G having universal cover T, the universal cover T minimizes the (normalized) 
partition function \n.Pc{z) for all z > 0, in particular it minimizes the (normalized) number 
of matchings Y2k m k(G). Of course, T being infinite, the normalized partition function needs 
to be defined properly and this can be done thanks to the local weak convergence nn. Indeed 
in the proof of Theorems [3] and [0 we will prove: 



Proposition 2. Let G be a finite bipartite graph and T be its universal cover. Let Q be the set 
of finite bipartite graphs with universal cover T. We have for all z > 0, 

inf -j— In P G '{z) = —max <j ( ^ x e ] Inz + Sq(x.) > . (12) 

G'egv(G') v(G) xeM(G) 1 \~^ I I 

Moreover, the sequence ( G n ) defined in Theorem^ or the sequence L n {G) converge in the local 

weak sense to T and achieves the bound (QU) in the limit when n tends to infinity. 


Note that G\, G 2 € Q if and only if G\ and G 2 have a common finite cover which is a result 
proved in (26]. Also, the right-hand term in (1121) is an invariant of Q\ this expression will be the 
same for any graph belonging to Q. Indeed our proof will proceed by computing its value thanks 
to the following message passing algorithm: to each edge uv € E and time step t, we associate 
two messages y t u ^. v (z) and obtained by setting y®-> v = Vv-m = 0 and for t > 0, 


„*+i _ 


(*) = 


^ "f" 5 ~2wGdu\v Vw^u( z ) 


(13) 


We show that as t —>• 00, these iterations will converge to a limit y u ^. v (z) and that x uv (z) = 
yu^ v (z)yv-,u{z) so ^ veg maximization in (IT!?h (see Propositions [SJ E] [7]). Since the recursion 
(1131) is local, the messages obtained after t iterations are the same as the one computed on the 
computation tree of the graph at depth t. We are able to show that these recursions on infinite 
trees still have a unique fixed point (see Theorem [9]) so that this fixed point should be the limit 
(as the number of iterations tend to infinity) obtained by our algorithm runned on the original 
graph G. Note that the recursion (1131) is well-known and first appeared in the analysis of the 
monomer dimer problem [233. H is a iso used to define a deterministic approximation algorithm 
for counting matchings in [6]. Indeed, their analysis directly implies that the convergence of our 
algorithm is exponentially fast in the number of iterations t. Note however that the recursion 
used in [6] corresponds to messages sent on the tree of self-avoiding paths. Instead, we use the 
tree of non-backtracking paths. The tree of self-avoiding paths is finite and depends on the root 
whereas our tree is the universal cover of the graph. Also, [6] directly implies the convergence 
of our message passing algorithm, it does not give any indication about the value of the limit. 


At this stage, we should recall that computing the number of matchings falls into the class 
of ^P-complete problems as well as the problem of counting the number of perfect matchings 
in a given bipartite graph, i.e. computing the permanent of an arbitrary 0 — 1 matrix. By 
previous discussion, we see that if the graph is locally tree like, then the tree of self-avoiding 
paths and the universal cover are locally the same, and one can believe that our algorithm will 
compute a good approximation for counting matchings. This idea was formalized in [43j and 
proved rigorously in im for random graphs. Our Theorem [5] shows that these results extend to 
random lifts. The lower bound in ([2]) is called the (logarithm of the) Bethe permanent in the 
physics literature @21 nans]. Similar ideas using lifts or covers of graphs have appeared in the 
literature about message passing algorithms, see [Ml M] and references therein. We refer to [25] 
for more results connecting Belief Propagation with our setting. 


3.2 Matching measure and spectral measure of trees 

We now relate our results to the matching measure used by Csikvari in nn and show how 
our results allow us to compute spectral measure of infinite trees. The matching polynomial is 
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defined as: 


V (G) 

Qg(z ) = £(-l ) k m k (G)z n ~ 2k = z n P G (-z~ 2 ). 
k=0 

We define the matching measure of G denoted by pc as the uniform distribution over the roots 
of the matching polynomial of G : 


PG 


1 

HG) 


u(G) 

^2 Sz i’ 


i=l 


(14) 


where the zfi s are the roots of Qg • Note that Qg(—z) = (—1 ) n Qc{z) so that pc is symmetric. 
The fundamental theorem for the matching polynomial is the following. 

Theorem 6. (Heilmann Lieb 123$ ) The roots of the matching polynomial Qg(z) are rea l an d in 
the interval [—2 \/Dg — 1, 2^/Dq — 1], where Dq is the maximal degree in G. 


In particular, the matching measure of G is a probability measure on M. Of course, the 
polynomials Pg(z) or Qg(z) contains the same information as the matching measure pc- We 
can express the quantity of interest in term of pg (see Lemma 8.5 in [23], [TJ or mr- for z > 0 , 


1 zP' G {z ) 
v{G) P G (z) 

HG) 

v(G) 


\ / TT^ dpa{x) ’ 

1(i-pg«o})). 


As explained above, Csikvari m uses this representation and the fact that for a sequence 
of d-regular graphs converging to a d-regular tree, the limiting matching measure is given by 
the Kesten-MacKay measure, to get an explicit formula for the limiting partition function. Our 
approach relies on local recursions instead of the connection with the matching measure. Since 
we are able to solve these recursions, we get the following result for the limiting matching 
measure. 

Theorem 7. Let G be a finite graph and T{G) be its universal cover. For any sequence of 
graphs ( Gi) with maximal degree Dg and with local weak limit T(G), the matching measure pGi 
of the graph Gi is weakly convergent to some measure Pt(G) defined by the formula for z > 0, 

/ ^ (l>w) • 

where the vector x(z) is the unique maximizer of ^^(x, z) in FM(G). Moreover, we have 
Pt(G )({°}) = 1 - 2 i^rr 

Note that our theorem gives a generating function of the moments of Pt(G) since for z > 0 
sufficiently small, we have: 

r -,\2 °° r 

J 1 I z \‘>df J ' T (G)W = y~!(~ 1 )' +1 ^ J ^ 2l dpT{G)W, 

Z— 1 
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and the series is convergent since the support of all the p Gi and hence of Pt(g) is contained 
in [—‘1\JDq — 1, 2 \JDq — 1], As shown by Godsil in mi, for finite trees, the spectral measure 
and the matching measure coincide, this is still true for infinite trees mmm- In particular, 
the moments of the matching measure f A 2t dpj-( G )(A) can be interpreted as the average number 
of closed walks on T(G) where the average is taken over the starting point of the walk (see 
Proposition [TT] for a precise definition of the random root as the starting point of the walk). 

To be more precise, for a finite graph G, we denote by Ai < ■ ■ ■ < A v (g) the real eigenvalues 
of its adjacency matrix and we define the empirical spectral measure of the graph G as the 
probability measure on R: 

x v(G) 
v 7 i=l 


The following theorem follows from mmm (see also Chapter 2 in [29] ) 

Theorem 8. Let G be a finite graph and T(G) be its universal cover. For any sequence of 
graphs (Gfi with maximal degree Dq and with local weak limit T{G), the spectral measure p Gi 
of the graph Gi is weakly convergent to the measure Pr(G) defined in Theorem [7| Moreover for 
all x € R, we have lim^oo p Gi ({x}) = p T (G)({ x })- 


In particular, we can apply this theorem to characterize the limiting spectral measure of 
random lifts L n {G ) as a function of the original graph G. For the atom at zero, we have: 


lim MLn(G )({°}) = 1-2 


v*(G) 

v(G) ' 


and Theorem [7] allows us to get the moments of the limiting measure of Pr n (G) ■ 


(15) 


4 Proofs 


4.1 Statistical physics 


To ease the notation, we will consider a setting with a weighted graph G = (V, E ) with positive 
weights on edges {9 e }e&E- Taking a bipartite graph G and 9 e = 1 for all e E E, we recover 
the framework of Section O To recover the more general framework of Section 12.21 consider 
the bipartite graph described by the support of A seen as an incidence matrix and for each 
e = (fii) € E, define 9 e = a % ^. 

We introduce the family of probability distributions on the set of matchings in G parametrised 
by a parameter z > 0: 


Tg( b) 


Be 


n 


eGB 1 


Pc(z) 


(16) 


where P G (z) = £b B ' : ILsb YlveV 1 (Eeeft; B e<^)= EfcS w k(G)z k , with 


u ’fc( G ) = e n 

(B:E e B e =k} eeB 
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where the sum is over matchings of size k. Note that we have w k (G) = per fc (A). Note also that 
when z tends to infinity, the measure /Iq converges to the measure: 


A#(B) 


Ilee E 

P er t/(G)(^) ' 


which is simply the uniform measure on maximum matchings when 9 e = 1 for all edges. In 
statistical physics, this model is known as the monomer-dimer model and its analysis goes back 
to the work of Heilmann and Lieb [23]. 

We define the following functions: 

U S G {Z) = -X>G(£e = 1), 

e&E 

U e G (z) = 5>fc(S e = l)ln(9 e , 

e&E 

S G (z) = -5>&(B)ln^(B). 

B 

Note that when 9 e = 1, we have Uq(z) = 0 and U G is called the internal energy while S G is the 
canonical entropy. We now define the partition function <1 >g(a) by 


$g0) = -Uq{z) In z + Uq(z) + S G (z). 


A more conventional notation in the statistical physics literature corresponds to an inverse 
temperature /3 = In z. Note that with our definitions, the internal energy Uq(z) is negative, 
equals to minus the average size of a matching sampled from fi z G . This convention is consistent 
with standard models in statistical physics where the low temperature regime minimizes the 
internal energy, i.e. in our context maximizes the size of the matching. A simple computation 
shows that: 

®g{z) = In P G (z) and, & G (z) = 

Lemma 1. The function U G (z) is strictly decreasing and mapping [0, oo) to (— v(G),0\. 


Proof. We have — U G (z ) = ^ w k(G)z k /P G (z) so that taking the derivative and multiplying 

by z, we get: 


<m\*) 


'Ek k2 ' W k( G )z k fJ2k kw k( G )z‘ 


Pg(z ) 


E *- 


Pea) 

J2e?we(G)z e \ 2 w k (G)z 


Pg(z) 


Pg{z) 


> 0 . 


□ 

We define r = t(G) = 2 v(G)/v(G) which is the maximum fraction of nodes covered by a 
matching in G. Note that t(G') < 1 and t(G) = 1 if and only if the graph G has a perfect 
matching. For t e [0, r), we define zt(G) € [0,oo) as the unique root to U G {zt{G)) = —tv(G)/2 . 
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Note that t i->- zfiG) is an increasing function which maps [0, r) to [0, oo). The function £ G (t) 
is then defined for t E [0, r) by: 


S G (i) 


S G (z t (G)) + U e G (z t (G)) 
v(G) 


(17) 


and £ G (t) = — oo for t > r. 

Proposition 3. Fort < r, u>e have £g(t) = — ^ln^(G). The limit limt^. T X G (t) exists and we 
define £ g (t) = lim t _> T S G (t) = In w v ^(G). 


Proof. We have for t < r, £ G (t) = ;^yy In P G (zt) — t/2\nzt, so that taking the derivative with 
respect to t, we get: 


s' G (t) 


zi _ p h( z t) \ _ 

2z t v(G)P G (zt)) 2 


Since Uq(z) = and C7 G (^t) = —tv(G)/ 2 , we get X' G (f) = —^lnz t . For t large enough, 

we have ^ > 1 and the proposition follows. □ 


The following proposition is proved in [141 for unweighted graphs (see Proposition 2.1(g)) 
but the proof is the same in the weighted case. We include it her for convenience. 


Proposition 4. If for some graphs G\ and G 2 , we have for every z >0, 

$Gi(z) $Ch(z) 

v{G\) ~ v(G 2 ) ’ 


then 


S Gl (t) > S G2 (i) 


for all 0 < t < 1. 

Proof. The assumption ensures that . Moreover if y = j , then 

ln^GiiCGi) ^ In w u ( G2 )(G 2 ) 
viGJ ~ v(G 2 ) ' 

Hence the statement is trivial for t > 2v(G 2 )/v(G 2 ). We consider now t E [0, 2v(G 2 )/v(G 2 )). 
Note that S Gl (0) = £ G 2 (0) = 0. The derivative of £ Gl (t) — S Ga (t) for t < 2u(G 2 )/v(G 2 ) is 

(ln^(Gi) - In z t (G 2 )) 

Assume this derivative is 0 at to, then we have zt 0 (G 1 ) = zt 0 {G\) = zo and then 
Sgi(~o) lnP Gl ( 2 0 ) *o , . lnP G2 (^ 0 ) to, S G2 (z 0 ) 

~Fov “ “W " 2 ln; 20 - “W ~ 2 ln 

Hence the minimums of S Gl (t) — S Ga (t) on [0,2n(G 2 )/v(G 2 )) are non-negative. □ 
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4.2 Local recursions on finite graphs and infinite trees 


Let G = ( V., E ) be a (possibly infinite) graph with bounded degree and weights on edges {0 e } e& E- 
We introduce the set ^ of directed edges of G comprising two directed edges u —» v and v —» u 
for each undirected edge (uv) € E. For if € 2^, we denote by —if the edge with opposite 
direction. With a slight abuse of notation, we denote by dv the set of incident edges to v E V 
directed towards v. We also denote by dv\u the set of neighbors of v from which we removed u. 
We also use this notation to denote the set of incident edges to v directed towards v from which 
we removed u —» v. 

Given G, we define the map TZg : (0,oo)^ —> (0,oo)^ by TZg(sl) = b with 


b, t — w — 


1 + Ylw£du\v @wu a w^u 


with the convention that the sum over the empty set equals zero. We also denote by TZ U ^ V : 
(0,oo) 9n \' !) —» (0, oo) the local mapping defined by: b u ^ v = lZ u -> v (a) (note that only the co¬ 
ordinates of a in du\v are taken as input of TZ U -^ V ). Comparisons between vectors are always 
componentwise. 


Proposition 5. Let G be a finite graph. For any z > 0, the fixed point equation y(z) = 
zTZc{y{z)) has a unique attractive solution y(z) € (0,+oo)^. The function z i->- y (z) is in¬ 
creasing and the function z ^ is decreasing for z > 0. 


Note that the mapping zTZg defined in this proposition is simply the mapping multiplying 
by z each component of the output of the mapping TZq (making the notation consistent). 


Proof. This result is proved for the case 6 e = 1 for all edges in [2SJ (see also [.‘Fij i and the proof 
extends to this setting. □ 


We define for all v € V, the following function of the vector (y-p, 


if € dv), 


Vv(y) 


£ 


9ey-?'R'--t( y) 

1 + 0 e y-^TZ_^{y) 


1 


(18) 

(19) 


Clearly from (1191) , we see that T> v is an increasing function of y and the proposition below follows 
directly from the monotonicity of y(z) proved in Proposition [5j 


Proposition 6. Let G = (V,E) be a finite graph and y(z) be the solution to y(z) = z7ZG(y(z)). 
For any v E V, the mapping z i->- T> v (y(z)) is increasing and V v (y(z)) = J2e£dv x e( z )’ where 


x e (z) 


e e y^{z)y_^{z) 
z + 8 e yi?(z)y-i?{z) 


( 20 ) 


We denote by x(z) = (x e (z), e € E) the vector defined by F2U\). then x(z) € FM(G) and we 
have: 


lim y^ V v (y(z)) = 2v*{G). (21) 

vev 
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Proof. The only non-trivial statement in the above proposition is the value of the limit in (1211) . 
In the case 9 e = 1, it follows from Theorem 1 in [28] and the proof carries over to the case 
9 e > 0 . □ 


For a finite graph G = (V,E) with weights on edges {9 e } e& E, we define for x £ FM(G) 
defined by (J5]) and z > 0, 

UB(X) = -^Xe, 

eOE 

Sg( x ) = y> e ln^ + (l- x e ) ln(l -x e ) - | 1 - Yl x e j In ( 1 - ^ x e J , 

e£E e vGV V e£dv J V e£dv / 

$g( x ,-) = -U§(x)]nz + S|(x). 

We denote by x(z) the vector dehned by (1201) in Proposition [ 6 ] where y(z) = zlZc T (y{z)). Note 
that 

D|(xW) = P- E V„(y(z)), ( 22 ) 

V&V 

so that by Proposition [ 6 ] the mapping z i->- Uq{x(z)) is decreasing from [0,oo) to (—i/*(G),0]. 
Thus, we can define z^ as the unique solution in [0, oo) to 


U8{x(z?)) 


B\ 


Similarly as in (ED, we define 


tv(G ) 
2 


for t < t*(G) 


2 v*[G) 

W 


E|(0 = for * < r*(G). 

Note that we have t*(G) > r(G) with equality if G is bipartite. 

Proposition 7. Recall that x(z) £ FM(G ) is defined by l2R) . We have for any z > 0, 

sup $ G (x;.z) = ^g(x(^);z), 
xeFM(G) 


and for t < t*{G), 


£g(t) = 


max 


W(G) xGFM t „ (G)/2 (G) 


5^ 


(x), 


where FM t is defined in [W and where the maximum taken over an empty set is equal to —oo. 


Proof. The first statement is proved in [28] for the case where 9 e = 1 but extends easily to 
the current framework. For the second statement, note that for any x £ FM tv iQy 2 {G) with 
t < t*(G), we have 


$g( x 




tv(G ) 


In m + Sg(x) < = 






tv(G) 


Inz? + S%(x(zf)) 


B\ 


By definition, we have x(zf) £ M to(G)/ 2 (G ! ), so that max xeM ^ (G)/2 ( G) S^(x) = S^fx(zf)). □ 
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We now extend Proposition [5] to infinite trees: 


Theorem 9. Let T = (V,E) be a (possibly infinite) tree with bounded degree. For each z > 0, 
there exists a unique solution in (0, oo)^ to the fixed point equation y(z) = zTZriyiz)), i.e. such 
that 


y u -+v{z) 


z 


i + E 


w£du\v 


^wuHw—^uy^) 


(23) 


Proof. First note that any non-negative solution must satisfy y u ^ v (z) < z for all (uv) E E. The 
compactness of [ 0 , z\^ (as a countable product of compact spaces) guarantees the existence of 
a solution by Schauder fixed point theorem. 

To prove the uniqueness, we follow the approach in [ 6 ]. First, we define the change of 
variable: h u ^ v = — In Vu ^A z ) so that (1231) becomes: 


= In | 1 + z ^2 6 wu e 

w£du\v 


Hy_ 


(24) 


We define the function / : [0, +oo) d e-)- [0, oo) as: 

/ k 


/(h) = ln 1 + zJ] 




where the parameters k, hi, 9 t . 9) and z are fixed and d = fF k j= \ ki. 

Iterating the recursion (1241) . we can rewrite it using such a function / so that uniqueness 
would be implied if we show that / is contracting. 

For any h and h 7 , we apply the mean value theorem to the function /(ah + (1 — a)h 7 ) so 
that there exists a E [0,1] such that for h Q = ah + (1 — a)h 7 , 


|/(h) - /(h')| = |V/(h a )(h - h')| < ||V/(h Q )|| Ll ||h - h 7 


A simple computation shows that: 


l|V/(h)|| Ll 


g. 

z 2^i=i u ', _ h i 

l+zT,7=x6}e > 


i + *EiU 


Let Aj = ^1 + z Ej=i Q) e hj ^j j then we get 

3 E£= 1 e ii A i " A i) = x _ 1 + g ELl Mi 

1 + z Ei=1 9'iAi l + Z Ef=l @i-Ai 

By taking the partial derivatives, we note that this last expression is maximized when all Aj are 
equal. Then the solution for the optimal Aj reduces to a quadratic equation with solution in 
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[0, +oo) equals to A* = ^ 1+ z Z q 1 , where 0 = Yli=i @i- Substituting for the maximum value, we 
get for any real vector h, 


l|V/(h)|| Ll < 1 


2 

y/TTzQ+l' 


□ 


4.3 2-lifts 

If G is a graph and v G V(G), the 1-neighbourhood of v is the subgraph consisting of all edges 
incident upon v. A graph homomorphism n : G' —> G is a covering map if for each v' G V(G'), it 
gives a bijection of the edges of the 1 -neighbourhood of v' with those of v = G' is a cover 

or a lift of G. If edges of G = (V. E) have weights 9 e then the edges of G' = (1 7',E') will also 
have weights with 9 e i = 9 n r e i\. Note that the definition of 2-lift for matrices given in Section T2.21 
is consistent with the definition of 2-lift for graphs by identifying the matrix A as the weighted 
incidence matrix of the bipartite graph. 

Proposition 8. Let G be a bipartite graph and H be a 2-lift of G. Then Pg(z) 2 > Ph(z) f or 
z > 0, So(t) > £#(i) for t G [0,1] and v{H) = 2 v{G). 

Proof. The proof follows from an argument of Csikvari nu. Note that GUG is a particular 2 -lift 
of G with Pgug{z) = Pg{z) 2 . To prove the first statement of the proposition, we need to show 
that for any 2-lift H of G, we have: Wk(G\JG) > Wk(H). Consider the projection of a matching 
of a 2 -lift of G to G. It will consist of disjoint union of cycles of even lengths (since G is bipartite), 
paths and double-edges when two edges project to the same edge. For such a projection R = 
R\ U i ?2 C E where i ?2 is the set of double edges, its weight is H[ egRi ^ n e gi? 2 ®e- Now for such 
a projection, we count the number of possible matchings in GllG: n^(G U G) = 2 k ^ R \ where 
k(R) is the number of connected components of R\. The number of possible matchings in H is 
ur(H) < 2 k ( R > since in each component if the inverse image of one edge is fixed then the inverse 
images of all other edges is also determined. There is no equality as in general not every cycle 
can be obtained as a projection of a matching of a 2-lift. For example, if one considers a 8 -cycle 
as a 2-lift of a 4-cycle, then no matching will project on the whole 4-cycle. 

Hence we proved that Wk(GL)G) > so that Pg(z) 2 > Ph(z ) for z > 0 and the second 

statement follows from Proposition [4j For the last statement, since Pg(z) 2 > Ph{z), we have 
2 v{G) > n(H) but the opposite inequality is true for any graph G since a maximum matching 
in G can be lifted to a matching in H with size twice the size of the original matching. □ 

Given a graph G with a distinguished vertex v G V, we construct the (infinite) rooted tree 
(T(G),v) of non-backtracking walks at v as follows: its vertices correspond to the finite non¬ 
backtracking walks in G starting in v, and we connect two walks if one of them is a one-step 
extension of the other. With a slight abuse of notation, we denote by v the root of the tree 
of non-backtracking walks started at v. Note that also we constructed T{G) from a particular 
vertex v, this choice is irrelevant. It is easy to see that T(G) is a cover of G , indeed it is the 
(unique up to isomorphism) cover of G that is also a cover of every other cover of G. T(G) is 
called the universal cover of G. 
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Since the local recursions are the same for both TIt(G) and Hg and since there is a unique 
fixed point for both zTZ t ^ and zTZq, the proposition below follows: 

Proposition 9. Let G be a finite graph and T(G ) be its universal cover and associated cover 
7 r : T(G) —>• G. By Propositions and\5i we can define: 

y{z) = zU T (G)(y(z)) , and, y(z) = zU G (y(z)). 

We have vr(y (z)) = y(z), i.e. y-f{z) = y^){z). 

4.4 The framework of local weak convergence 

This section gives a brief account of the framework of local weak convergence. For more details, 
we refer to the surveys mm- 


Rooted graphs. A rooted graph (G, o ) is a graph G = (V, E) together with a distinguished 
vertex o G V, called the root. We let G * denote the set of all locally finite connected rooted 
graphs considered up to rooted isomorphism, i.e. (G,o) = ( G',o') if there exists a bijection 
7 : V —>• V' that preserves roots ( 7 ( 0 ) = o') and adjacency ({i,j} G E ( 7 ( 2 ), 7 (j)} G E'). 
We write [G, o\h for the (finite) rooted subgraph induced by the vertices lying at graph-distance 
at most h G N from o. The distance 

dist ((G, o), (G', o’)) := —j— where r = sup {h G N: [G, o\ h = [G 1 , o’] h } , 

' 1 + r 

turns G* into a complete separable metric space, see [ 2 ], 

With a slight abuse of notation, {G, o) will denote an equivalence class of rooted graph also 
called unlabeled rooted graph in graph theory terminology. Note that if two rooted graphs are 
isomorphic, then their rooted trees of non-backtracking walks are also isomorphic. It thus makes 
sense to define (T(G),o) for elements (G,o) G £?*. 

Proposition 10. For any graph G = (V,E), there exists a graph sequence {G n } nS N such that 
Gq = G, G n is a 2-lift of G n -\ for n > 1. Hence G n is a 2 n -lift of G and we denote by n n : G n —>• 
G the corresponding covering. For any uGf, if v n G 7 r“ 1 (u), we have (G n ,v n ) —» ( T{G),v ) in 

G *. 

Proof. The proof follows from an argument of Nathan Linial [30] . see also [14J. 

A random 2-lift H of a base graph G is the random graph obtained by choosing between the 
two pairs of edges ((u, 0), (v, 0)) and ((u, 1), (v, 1)) G E(H) or ((u, 0), (v, 1)) and ((u, 1), (v, 0)) G 
E(H ) with probability 1/2 and each choice being made independently. 

Let G be a graph with girth 7 and let k be the number of cycles in G with size 7 . Let X 
be the number of 7 -cycles in H a random 2-lift of G. The girth of H must be at least 7 and 
a 7 -cycle in H must be a lift of a 7 -cycle in G. A 7 -cycle in G yields: a 27 -cycle in H with 
probability 1/2; or two 7 -cycles in H with probability 1/2. Hence we have E[A] = k. But GuG 
(the trivial lift) has 2k 7 -cycles. Hence there exists a 2-lift with strictly less than k 7 -cycles. By 
iterating this step, we see that there exists a sequence {G n } of 2-lifts such that for any 7 , there 
exists a n( 7 ) such that for j > n( 7 ), the graph Gj has no cycle of length at most 7 . This implies 
that for any v G V and Vj G we have DIST ( (Gj , v 3 ) , (T(G),v)) < ^ and the proposition 

follows. □ 
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Local weak limits. Let V(Q+) denote the set of Borel probability measures on £/*, equipped 
with the usual topology of weak convergence (see e.g. 0). Given a finite graph G = (I/, E), we 
construct a random element of G* by choosing uniformly at random a vertex o E V to be the 
root, and restricting G to the connected component of o. The resulting law is denoted by U{G). 
If {G n } n >i is a sequence of finite graphs such that {U{G n )}n>\ admits a weak limit £ E V{Q *), 
we call £ the local weak limit of {G n } n >i. If ( G , o) denotes a random element of £/* with law £, 
we shall use the following slightly abusive notation : G n (G, o) and for / : £/* —> M: 

E ( g,o)[/(G,o)] = f f(G,o)dC(G,o). 


As a direct consequence of Proposition [TOl we get: 

Proposition 11. For G = (V,E), let be the sequence of 2-lifts defined in Proposition 

m Then G n ( T(G),o ) where T(G ) is the universal cover of G with associated cover n : 
T{G ) — >• G and o is the inverse image of a uniform vertex v of G, o = 7r _1 (u). 

We now state the corresponding well-known result for random lifts: 

Proposition 12. For G = (V,E), let L n (G) be a random n-lift of G. Then L n (G) (T(G),o) 
a.s. where T(G) is the universal cover of G with associated cover 7r : T(G) —» G and o is the 
inverse image of a uniform vertex v of G, o = 7r _1 (u). 


We are now ready to use the results of the above sections. The existence of the limits for the 
partition function, the internal energy of the monomer-dimer model is known to be continuous 
for the local weak convergence (in a much more general setting than here) [2f> . II, 271 fT] but the 
explicit expressions given in the right-hand side below are new. 


Theorem 10. Let G be a finite graph and T{G) be its universal cover. Let (G n ) n >i be a sequence 
such that G n (T(G),o). We denote by x(z) the vector defined by EM) in Proposition^ where 
y (z) = zTZc(y{z)). Then we have as n —» oo, for z > 0, 


lim 

n—>• oo 


lim 

n—»• oo 


lim 

71—> OO 


1 

Vn 


lim 

n^-oo 


1 


AGr 


| V n 
I In Pc n (z) 
1 r uzjz) 


Vn 



lim T Gn {t) 
n—> oo 


v(G) 

1 

Ag) 

i 

Ac) 

i 


v*(G), 


®g(Az),z) 


U g A(z)), 

SgA(z)), 


v(G) 

Sg(*), f or t < t *(G). 


(25) 

(26) 

(27) 

(28) 
(29) 


Proof. In [231 ED, it is shown that the root exposure probability satisfies (with our notation): 

Tu —— 


1 Z ^2w£du\v ^wuTw^u{z) 
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Hence we can use directly results from HD by the simple change of variable: y u ^. v (z) = zr u ^. v (z). 
In particular Theorem 6 in HD implies that 


2 E (T(G),o) 

2 E (T(G),o) 

^ e (t(g),o) \Po(y(z ))\, 

and (|27|) follows from Propositions [9] and [6j (1251) follows by taking the limit z —>• oo as shown in 
Theorem 11 in HD and (12T1) in Proposition [6] 

We now prove (1261) . We start by noting that $'q(z) = Ug ^ so that the convergence of 
t^-t In Pq u (z) follows from (1271) and Lebesgue dominated convergence theorem (see Corollary 7 
in HD)- We only need to check the validity of the right-hand side expression in (1261) . 

Note that, we have with 0 m i n = min e 0 e > 0 and 0 max = nrax e 6 e > 0 


1 

1 + z Ei tGdo°er^(z)_ 
Y,^ & do e ey^{z) 

- 1 + ^eUt{z)_ 


1 

TT77 U G„( Z ) = 

n—>• oo | v n | 


1 In PGn{z) > _ JJS / X l-E'nl ln^min 


1141 lnz <jrnV ' ’ 114| lnz ’ 

and since the number of matching is upper bounded by 2^ En ^, we have 


1 In Pg„(z) < s . . \E n \ ln0 max 
— G n ' / ' 


I E n | In 2 

114 | In 2 — u "' ' ' 1141 In z ' 1141 In 

Hence, taking first the limit n —>• oo and then the limit z —>• oo, we have 

1 In P Gn (z) v'{G) 


lim lim 

z—>-oo n—>-oo \ v n 


In z 


v(G) 


Since $>q(x.(z), z) ~ ^|^ylnz by Proposition [6] (note that Sq(x) is bounded), we only need 

U B (x(z)) 

to check that the derivative with respect to z of the right-hand term in (1261) is G , ~ . 
Lemma 2. In the setting of Proposition^ we have 

x e (z)(i -*.M) =>c xA 7 (30) 

V e’£du ) V e'&dv ) 


Proof. Note that ^2f £ g v Xf(z) = V v (y(z)), so that we have by (fl9l) 


1 - J2 x f (*) 

f&dv 


= 1 - 


Elfg dv 9 ^( Z ) 

l + 12~tedv 9 e y-£(z) 


1 + e hyp(z) 
it £dv 
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We have for e = ( uv ) E E, 


x e (z) = 


^eUu—¥v (^) 


Vv- 


^eVu—^vi^) 


and using the fact that y(z) = zTZa(y(z)), we get 

*e(*) = T -M = deVu ^ v ^ ( ' 1 “ S X f^ 

1 “T Z^wGdv v wvVw^v\Z) \ j edv 

i ^ ^w€du\v QwuUw->u(z) _ Z f i ^ 

J- Xp\z) — “ ^~ 7 “ — 7 r I J- / v £/(£ 


i + E O wu y w ^ u (z) y u ^ v (z) 


f&du 


and the lemma follows. 


□ 


Note that for e = (uu), we have 


d ^G , 

—= ln z + In 

ox e 



q<$>b d<& B D 

In particular, we have -^ s -(x(z)) = 0 by Lemma [2] and then -^f-(z) = —Uq(x.(z))/z and (1261) 
follows. Moreover (1281) follows from (l26l) and (l27l) . 

We now prove (1291) . Assume that there exists an infinite sequence of indices n such that 
zt(G n ) > zf + e. We denote zi = zf and Z 2 = zj 3 + e. We have for those indices: 

-W\ U& ~ M - £ -]h uiAzt(a ' )) = 5 - 

Then by the first part of the proof, we have — ^g„(zi) ~ (x(zi)) = | and 

~]h U Gn( Z ^ — > — ^q^Uq (x(z 2 )) > | by the strict monotonicity of z Uq(x(z)). Hence 
we obtain a contradiction. We can do a similar argument for indices such that zt{G n ) < zj 3 — e, 
so that we proved that zt(G n ) —> zj 3 . Then (1291) follows from the continuity of the mappings 
z i->- y(z) and x i->- Sq(x). □ 

Proposition 13. The function Sq (x) is non-negative and concave on FM(G). 


Proof. From Theorem 20 in [39], we know that the function 

/i(x) = - y~] Xj In Xj + y^(l - Xj) ln(l - Xj) 

i i 

- ^ x i j ln ^ H x ^j + x ^j ln 

is non-negative and concave on A k = {x E K fc , Xi > 0,39^1] Xi < 1}. Hence the function 
fif(x) = - ^ Xj ln Xj + y^(l - Xj) ln(l - xf) - 2 1 - ^ Xj In 1 - ^ Xj\ 
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is concave and non-negative on A k since 


ff( x ) 



where H(p ) = — plnp — (1 — p) ln(l — p) is the entropy of a Bernoulli random variable and is 
concave in p. The proposition follows by decomposing the sum in Sjf (x) vertex by vertex. □ 


4.5 Proof of Theorem |4j 

Corollary 1. Let G be a bipartite graph, then for any z > 0, 

$g{z) = lnPclz) > max &£(x.:z) 
y K xe FM{G ) GK 


and for t < t*{G), we have 


S G (t) > 


max 


v(G) *eFM h(G)/J (0 



Proof. We consider the sequence of graphs defined in Theorem [101 By Proposition [51 the se¬ 
quence {-|^-|4 > G„(^)}neN is non-increasing in n and converges to ^4>g(x(z),z) by Theorem 
[TOl Hence the first statement follows from Proposition [7] 

The second statement of Proposition [8] implies that the sequence {S(j n (i)} ne N is non¬ 
increasing in n and converges to E^(t) by Theorem [TO] and the last statement follows from 
Proposition [7] □ 

The final step for the proof of Theorem 0] is now a standard application of probabilistic 
bounds on the coefficients of polynomials with only real zeros [33] . 

Let k < iy(G) = u, t = and 2 = zt(G) such that Uq(z) = —tv(G)/2 = —k. For i < v, 
we define 

Wi(G)z i 
Ul ~ Pg(z) • 

By the Heilmann-Lieb theorem [23], the polynomial A(x ) = Y/i=o aiX% ^ as on ly rea ^ zeros , he. 
(ao, • • •, a v ) is a Polya Frequency (PF) sequence. Note that A(l) = 1 = ]>T a*. By Proposition 1 
in [33], the sequence (ao,..., a u ) is the distribution of the number S of successes in v independent 
trials with probability pi of success on the i-th trial, where the roots of A{x) are given by 
— (1 — pi)/pi for i with p % > 0. Note that E[S] = Yli = — Ug( z ) = k. 

We can now use Hoeffding’s inequality see Theorem 5 in [24]: let S' be a random variable 
with probability distribution of the number of successes in v independent trials. Assume that 
E[S] = up G [6, c]. Then 


P(Se[M)>E 

i=b 


Qp'd-rt 
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Hence, we have in our setting with b = c = k and p = L -: 


O'k — 
w’fc(G) > 

> 


^P k (i-pY (x - p) 

b v ,k(p) exp (v(G)T, G (t)) 

h,k(k/v) exp max S G (x) , 

\y^&FM tv ^ G y 2 (G) J 


where the last inequality follows from Corollary [0 

The case k = u is easy. Take t = 2u ^q^ with e > 0 and z = zt(G ) so that Uq(z) = 
—tv(G)/ 2 = —i/(l — e). We define the sequence of s as above. We now have EfS 1 ] = u(\ — e). 
We then have EfS 1 ] = Yi < va v + (1 — a v ){y — 1) = a u + v — 1, so that a v > 1 — ve and 


Wv(G){G) > (1 


ue) exp (v(G)Tia(t)) > (1 


ve) exp 


max 

6™in(G)/ 2 (G) 



Letting e —> 0 concludes the proof. 


4.6 Proof of Theorem [5] 

We start with a definition: the perfect matching corresponding to the edge e in G is called the 
fibre corresponding to e, which we denote by F e . 

We denote u n = v(L n (G)). If G is bipartite, we have by Theorem [3] for all k < u n , 


\nm k (L n (G)) > ( £) + xeM “£ (G)) S Ll.O)W- 


It is easy to see that 


max Sr /r>i(x) > n max Sr(x), 
xe M k (L n {G)) Ln ^ )K xeM fe/n (G) 


since to any x E M k / n . we can associate y € M k (L n (G)) by taking y e r = x e for all e' E F e 
Hence, we get 


1 i ( k 

- In m k {L n (G)) > max S G (x) + - In b v k ( — 

Tl x£Mj./ n (G) TL \F n 


Taking k = v n and letting n —>• oo, we get 


1 


lim inf -In m Un (L n (G)) > max S G (x) 
n ->oo n x£FM„. (G )(G) 


For the upper bound, we do not need to assume that G is bipartite as we have for all z > 0, 


- hr m Vn {L n (G)) < - In P Lri{G) (z) - — \nz. 
n n v ’ n 


23 



Hence letting n —>• oo, we get 


lim sup — In m Vn (L n (G)) < ®q(x.(z), z) - v*{G) \nz 

n—too IT, 


= sup 
xeFM(G) 


M? 


x e - iy*(G) + Sq (x) 


Taking now z —>• oo and noting that Yl e x e ~ v*{G) < 0, we have: 

lim sup — In m Vn (L n (G)) < sup Sq (x). 
n->oo n x.&PM(G) 
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